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SUMM.\RY 

Wh«n  a.  rocket  is  launched,   there  is  s  short  initial  period  of  ac- 
celeration during  which  the  rocket  is  unstable.    As  the  flight  velocity 
Increases,  the  aerodynaxnic  forces  acting  on  the  fins  and  stabilisers 
become  large  enough  to  give  stability.     Various  methods  have  been  ezn* 
ployed  to  stabiUae  the  rocket  during  this  launching  period.    Cxiide  radls, 
"aero  length"  launchers,  booster  rockets  which  produce  high  initial  ac- 
celeration, and  auto -pilot  controlled  noasles  are  typical  devices  that 
have  been  used. 

This  is  an  investigation  of  the  requirenrients  of  a  nozcle  control 
which  would  stabilise  the  rocket  during  the  launching  period.     The  con- 
figuration investigated  is  unique  in  that  the  nosale  of  the  rocket  is 
n>ounted  as  a  compound  pendulum,   and  the  movement  of  the  pendulum 
is  utilised  to  famish  the  signal  for  the  nossle  control  servo-nnechan- 
ism,   thereby  eliminating  the  need  for  gyroscopic  elements  in  the  con- 
trol system.     The  penduluzn  motion  of  the  nossle  caused  by  a  change  in 
flight  attitude  of  the  rocket  is  introduced  into  a  compxiter  which  pro- 
duces an  output  signal  proportional  to  the  attitude  of  the  rocket.    This 
attitude  signal  is  fed  back  to  the  nossle  control^  which  positions  the 
nossle. 

The  results  of  the  analysis  showed  that  the  rocket  was  unstable 
during  the  take-off  period  when  the  nossle  control  acted  on  the  rocket 
attitude  signal  alone .     Stability  over  a  narrow  range  of  feedback  gains 
was  indicated  for  the  system  using  a  simple  lead  circuit  as  a  nozzle 
control,  or  in  other  words,  when  the  nossle  control  acted  on  both  the 
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attitude  signal  and  the  rate  of  change  of  attitude  of  the  rocket.     The 
damping  characteristics  of  this  system  were  poor.     By  changing  the 
nosmle  control  function  to  include  a  response  to  the  acceleration  of  the 
rocket  attitude,   the  damping  characteristics  were  in\proved  and  the 
range  of  feedback  gains  was  widened. 
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I.     INTRODUCTION 

The   rocket   Btabllity  problem  has   existed  for  xnany  centuries. 
Historians  record  that  rockets  were  first  used  in  battle  by  the  Chinese 
in  1232  A.D.    3y  1400,   rockets  were  in  widespread  use  as  major  wea- 
pons by  the  warring  nations  of  Europe.    The  poor  stability  of  the  rocket 
in  flight  was   responsible  for   its   gradual   replacement  by  the   cannon 
which  gave  nu>re  accurate  trajectories.     In  the  late  1700's,   India  em- 
ployed rockets  with  improved  stability  characteristics   against  the 
British.    V/hen  I^ritish  soldiers  attempted  to  defeat  the  forces  of  Hydar 
Ali,  Prince  of  Mysore,  a  trained  corps  of  rocket  boxnbardiers  launched 
rocket  missiles  at  the  charging  British  cavalry,  inflicting  many  casual* 
tics.     The  Indian  rockets  were  made  of  iron  cylinders  about  eight  inches 
long  and  two  inches  in  diameter,  with  a  bamboo  shaft  about  ten  feet  long 
attached  to  give  stability. 

i\ji  English  military  officer.    Colonel  v/iUiam  Congreve,  becan^e 
interested  in  the   Indian  advances   in  rocket  dcvelopnnent,    and  worked 
to  further  improve  the   rocket  for  British  use.     Congreve  devised  a 
launching  system  consisting  of  a  tube  mounted  at  an  angle   on  a  tripod 
base,   similar  to  the  mortar.    Long  shafts  were  attached  to  the  rockets 
to  keep  thein  flying  in  a  straight  line.     The  British  employed  rockets 
against  Napoleon  at  Poulogne,   and  in  other  bombardments  with  great 
success. 

During  the   V/ar  of  1812,    /onericans  were  first  introduced  to 
rocket  warfare  when  the  British  bombarded  Bladenburg,  Maryland,  and 
subsequently  Ft.  McHenry  at  Baltinrxorc.    An  American  by  the  name  of 
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Williazn  Hale  improved  the  stability  of  the  Congreve  rocket  by  replac- 
ing the  >vooden  shaft  with  three  curved  fins  on  the  rear  part  of  the 
rocket  body.     The  exhaust  gases  were  deflected  by  the  fins,   causing 
the  rocket  to  spin  in  flight.     These  rockets  were  used  in  the  war  with 
Mexico,   and  were  found  to  have  superior  trajectories  to  the  stick  sta- 
bilixed  types.     Again  the  rocket  declined  in  popularity  partly  because 
of  its  inferior  trajectory  compared  to  that  of  the  gun.     No  other  out- 
standing developnnents  in  rocket  stabilisation  were  made  following 
Hale's  work  until  early  in  this  century. 

With  the  advent  of  the  airplane,  and  an  increase  in  knowledge  in 
the  field  of  aerodynamics,  the  stability  problem  of  the  rocket  in  flight 
became  defined  nvore  clearly.     The  use  of  tail  fins,   or  tail  fins  with  a 
shroud  ring  as  on  the  Baaooka,   and  later,  the  use  of  controlled  airfoils 
solved  the  in-flight  stability  problem  to  a  great  extent. 

In  the  first  part  of  the  rocket's  flight  path  there  is  a  destabilising 
effect  from  the  change  in  mass  of  the  rocket  as  fuel  is  burned.     The 
aerodynamic  forces  must  counter  this  tendency  in  order  to  stabilise 
the  rocket.     During  the  launching  period,  the  aerodynamic  forces  are 
too  small  to  counteract  the  effect  of  the  change  in  nnass,  and  the  rocket 
is  unstable  or  at  best  neutrally  stable.     Consequently,    stability  during 
the  critical  take-off  period  must  be  attained  by  some  means  other  than 
the  use  of  aerodynamic  surfaces. 

Various  means  have  been  employed  to  stabilise  the  rocket  during 
the  launching  period.    Guide  rails  to  hold  the  rocket  to  a  straight  flight 
path  until  it  has  gained  sufficient  speed  to  become  aerodynamically 
stable  have  been  used  successfully.     The  length  of  the  guide  rail  has 
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been  reduced  to  a  mlnixnuzn  in  some  applications  by  firing  the  rocket 
with  a  very  high  initial  acceleration.     The  Cerxnan  V-2  rocket  was 
stabilised  during  the  launch  by  an  autonnatic  pilot  which  controlled 
four  carbon  vanes  in  the  Jet  stream  behind  the  noasle.     The  American 
"Viking"  and  the  German  "Enaian"  employed  a  controllable  notizle 
positioned  by  an  automatic  pilot. 

A  unique  means  of  launching  a  large  rocket  with  a  booster  rocket 
mounted  free  to  swing  as  a  connpound  pendulunn  was  investigated  by 
J.  C.  Norris  in  1951  following  a  proposal  by  H.  S.  Tsien  of  Cal  Tech. 
(Ref .   1).    For  the  system  investigated,  the  booster  rocket  was  nvsunted 
below  the  nudn  rocket  on  gimbals  off  the  center  of  gravity  of  the 
booster  as  shown  in  Fig.    1.     A  deviation  of  the  main  rocket  from  its 
course  would  cause  zr^tion  of  the  booster  relative  to  the  nrvain  rocket 
causing  a  change  in  the  direction  of  the  thrust  axis.     It  was  hoped  that 
a  practical  configuration  could  be  devised  to  give  launching  stability. 
The  results  of  the  Investigation  were  somewhat  discouraging  in  that 
for  a  typical  vehicle -booster  combination  for  which  stability  was  indi* 
cated,  the  required  length  of  the  tail  boom  was  too  great  to  be  prac- 
tical. 

Since  the  basic  idea  of  the  pendtdum  mounted  thrust  unit  appears 
sound,  the  possibility  of  utilising  the  pendulum  in  conjunction  with  a 
servo -control  is  suggested. 

In  this  thesis,   the  reqxxirements  of  such  a  servo-control  system 
which  would  stabilise  the  rocket  during  the  launching  period  are  inves- 
tigated.    The  rocket  has  a  controllable  noasle  mounted  as  a  compound 
pendulum;  conceivably,   the  motor  and  noxsle  of  a  liquid  propellant 
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rocket  could  be  moiinted  ai  a  unit  as  shown  in  Fig.  2.     The   motion  of 
the  pendtdum  with  respect  to  the  rocket  is  fed  to  a  computer  which  cal< 
culates  the  attitude  of  the  rocket.     The  attitude  signal  from  the   com- 
puter is  fed  to  a  nossle  control.     The  aim  of  this  investigation  is  to 
determine  the  requirements  of  the  nosale  control  such  that  it  will  sta- 
bilise the  rocket  during  the  launching  period. 

The  equations  of  nu>tion  of  the  configuration  were  first  derived, 
assuming  all  aerodynamic  forces  to  be   negligibly  small   during  the 
launching  period.     Motion  in  one  plane  only  was  considered.    External 
forces  considered  were  the  rocket  thrust,  assumed  to  be  constant;  the 
Jet  damping  force  arising  from  rotation  of  the  Jet  stream;  and  the  gra- 
vity force.     The  effect  of  the  change  in  mass  of  the  rocket  resulting 
from  the  burning  of  fuel  was  considered.     The  displacements  of  the 
rocket  and  nossle  about  the  initial  condition  were  assumed  to  be  small 
quantities  for  the  period  of  the  stability  investigation. 

Following  tiie  derivation,   the   equations   involving  the   pitching 
motion  of  the  rocket  and  the  motion  of  the  nossle  were  transformed  to 
simple  algebraic  equations  by  use  of  the  Laplace  Transform.     From 
the  transformed  equations,  it  was  possible  to   solve  for  the  transfer 
functions  of  all  of  the  components  of  the  control  system,   as  shown  in 
Fig.  4,   except  that  of  the  noazle  control.     The  form  of  the  nossle  con- 
trol transfer  function  was  estimated,   and  the  stability  of  the  resulting 
system  was  investigated  by  use  of  the  Nyquist  Stability  Criterion.    The 
Evans  Root  Locus  plot  was  used  to  find  the  damping  coefficients  of  the 
modes  of  oscillation  of  the  system.     Successive  estimates  of  the  form 
of  the  nossle   control  function  were   made   in  order   to  improve   the 
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•tabllity  of  the  system. 

It  was  found  that  for  a  nossle  control  function  ecnploying  feedback 
of  the  rocket  attitude  signal  only,   the  system  was  unstable.     For  a  con- 
trol function  eznploying  feedback  of  both  the  rocket  attitude  and  rate  of 
change  of  attitude,   it  was  determined  that  a  stable  range   existed  be- 
tween upper  and  lower  linnits  of  feedback  gain;  the  linnits  depending  on 
the  constants  of  the  control  function.     The  damping  of  the  system  was 
poor. 

The  nossle  control  function  was  revised  by  introducing  additional 
lead  circuits,   or  sensitivity  to  the  accelerative  change   of  attitude  of 
the  rocket.     The  damping  was  improved  by  this  change,   and  the  range 
of  gains  was  widened. 
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U.     AN-^J^YSIS 
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Fig.  3.    Block  DlAgrazn  of  a  Feedback  Control  System 


A  typical  feedback  control    system  which  might  be  employed  to 
guide  a  rocket  of  the  controllable  nozzle  type  is  shown  in  Fig.  3  above. 
The  attitude  of  the  rocket,    as  indicated  by  the   angle   Oy  which  the 
rocket  naakes  with  the  vertical,  is  measured  by  an  attitude  gyro.     The 
gyro  signal  is  zr.atched  with  the  desired  rocket  attitude  signal  and  the 
error,   or  difference,   is  fed  to  a  noaale  control.      The  nossle   control 
applies  a  moment  A)1   to  the  nosale  which  results  in  a  relative  angular 
displacement     Oy    Oj     of  the  nosale  with  respect  to  the  rocket.     The 
thrust  axis  rotates  with  the  noaale,    causing  a  change   in  the   rocket 
attitude . 

Now  consider  the  case  of  a  rocket  with  the   n\otor   and  nossle 
mounted  as  a  compound  pendulum  as  shown  in  Figs.  1  and  2.    With  any 
change  in  rocket  attitude,  the  pendulum  moves  in  a  calculable  nnanner . 
From  a  knowledge  of  the  system  dynamics  and.  a  time  history  of  the 
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penduluzn  motion,  the  attitude  of  the  rocket  may  be  computed.    Suppose, 
then,   that  the  angular  displacement  of  the  nosale  from  the  center  line 
of  the  rocket,    f  v  ~  ^5    *   ^^  measured  and  fed  as  an  electrical  or  me- 
chanical signal  into  a  computer  with  an  output  representing  the  attitude 
of  the  rocket   ^5     .     The  gyroscopic  element  is  thereby  eliminated 
from  the  circuit. 
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Fig.  4.     Block  Diagram  of  a  Feedback  Control  System 

Employing  a  Connpound  Pendulum  Noszle 


The  block  diagram  of  the  proposed  feedback  control  circuit  as 
shown  in  Fig.  4  above  may  be  interpreted  as  follows. 

The  desired  input  of  the  coordinate  ^       ,   or  angle  which  the 
rocket  makes  with  the  vertical,   is  matched  with  a  feedback  signal  rep- 
resenting the  actual  angle  at  any  given  time.     The  difference  between 
the  desired  angle  (in  this  case  it  is  desired  that    0^=0    for  vertical 
takeoff)  and  the  computed  angle  is  fed  into  the  nossle  control  naechan- 
ism  as  an  error  signal,    fj^.  -  9 3  1\ie.  error  signal  initiates   the 

application  of  a  moment  ^77  to  the  nozzle  gimbals .     The   moment 
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contributes  to  the  nozzle  dynamics  to  yield  a  certain  position  of  the 
nozzle  relative  to  the  rocket,   which  is  represented  ^Y  ^4  ~  %-x    •    ^^^ 
thrust  line  is  rotated  off  the  center  line  of  the  rocket  by  movement  of 
the  nozzle;  this  contributes  a  moment  and  side  force  to  the  overall 
rocket  dynamics,   resulting  in  a  final  rocket  position. 

The  functians     F2  .    P3  »   *^^    f^   appearing  in  the  block  diagram 
of  Fig.  4  are  transfer  functions  derived  from  the  equations  of  motion 
of  the  configuration  which  have  been  transformed  by  the  use  of  the 
Laplace  Transfornnation.     The  transfer  functions  are  linear  operators 
and  can  be  defined  by  inspection  of  the  block  diagram   of   the    system 
(i.e..  It-       Tj  .      rj"!^."     Q4-Q3       )•     The  no»»le  control  trans- 

fer function    p,    is  the  unknown  quantity  which  is  the  subject  of  this  in- 
vestigation. 

The  e.^uation  relating  the  output  of  the  system  to  the  input  n\ay 
be  derived  as  follows: 
Following  the  circuit  from    Q^  ^^    to      Q3  g^-p 

Around  the  loop 


Q3   =lo,,K-  Q,]f,FkF^ 


0     -       Ft  Fi  Ph  Q3  in 

Subatitutlng  for    Q-,  in  tlio  first  equation 


Q3  OUT  _   Output  _ 


F.  Ft  r4 
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The  System  Transfer  Function  therefore  becoznee: 

K     r     —  O^TFUT    _  Pi  l^t  '3 

^s^s  "  Input  [  +  p,  p^  P4 

The  ttability  analysis  thus  involves    1)  derivation  of  the  equations 
of  motion  for  the  rocket  and  movable  noaale  configuration,    2)  conver- 
sion of  these  equations  by  the  Laplace  Transform  to  solve  for  the  trans- 
fer functions    f^     «    f\    >   itnd    p^    ,  and  finally,    3)  synthesis  of  a  noaale 
control  function    p     ,   and  investigation  of  the  complete  system  trans- 
fer function    K.  G      ^^^  ^^*  stability  characteristics. 

The  Nyquist  Stability  Criterion  and  the  Evans  Root  Locus  Method 
are  applied  in  the  stability  investigation.     Refs.   (2)  and  (3)  contain  ex- 
planations of  the  use  of  these  and  other  methods  of  analysis  applicable 
in  the  field  of  stability  and  control. 
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m.     EQUATIONS  OF  MOTION 


The  assumed  rocket  configuration  with  all  external  forces  to  be 
considered  is  shown  in  Fig.   5.     .-.  schennatic  diagram  of  the  coordinate 
system  and  dimensions  of  the  configuration  are  shown  in  Fig.  6. 

Frictional  and  aerodymunic  forces  are  neglected.     The  change  of 
nnass  of  the  main  rocket  body  with  burning  of  fuel  is  considered.     Mo- 
tion  in  one  plane  only  is  considered. 

It  is  convenient  to  enr^ploy  coordinates  Q      through    Q     in  the 
derivation,   but  of  these,   only  four  coordinates.    Q       through    Q     .   are 
independent.     The  coordinates     ci     and     (^    measure  the  horizontal  and 
vertical  distances  respectively  to  the  center  of  gravity  of  the  entire 
configuration,   and  have  the  dimensions  of  length.     Coordinates    Q.    and 

Q     measure  the  angiilar  displacement  of  the  main  rocket  center  line 
and  the  noaale  line  respectively  from  the  vertical,  and  have  the  dimen- 
sions of  the  radian. 

The  equations  of  motion  in  the  four  degrees  of  freedoR\,  namely 
in  the  directions  of  coordinates  Q  ,  a  ,  q  ,  and  Q  are  derived 
using  the  generalised  L*agrangian  momentum  equation: 


d    /3T   \         dJ 


=  a.^  n» 


Elimination  of  Extraneous  Coordinates 

By  taking  moments  about  the  center  of  gravity  of  the  system  as 
shown  schematically  in  Fig.  7,  the  relation  between  the  system  center 
of  gravity  amd  the  centers  of  gravity  of  the  onain  rocket  body  and  the 
noszle  nnay  be  found: 
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Fig.   7.     Sketch  Showing  the  Relatloa  Between  the  Centers 
of  Mass,   Mj  mnd  M2,   and  the  Center  of  Gravity  of  the  System 


(2) 


From  the  geometry  of  the  eyetem  as  shown  in  Fig.  8  the  follow- 
ing relations  can  be  obtained: 
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Fig.   8.     Geometry  of  the  Configuration 


Equating  the  vertical  distances  between  nnass  centers: 
Equating  the  horizontal  distances  between  mass  centers: 
By  considering  sixnilar  triangles  in  Fig.   8   it  follows  that: 


(3) 


(4) 
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Substituting  £4.   (2)  into  (5),   the  result  la 

£41.   (2).   (3).   (4),   and  (6)  are  four  linear  expresaiona  containing 
the  four  extraneous  coordinates     ^^  .     ^    .     ^    ,     9  ^  *^n<l  the  principle 
coordinates     ^     ,     ^^  ,     <^     and     *^,  .     The  former  niay  be  expressed 
in  terms  of  the  principle  coordinates  as  follows: 
From  £q.   (3) 

From  Eq.   (4) 

Substituting  (8)  into  (2) 

Substituting  (7)  into  (6) 

Substituting  (9)  into  (8) 

Substituting  (10)  into  (7) 

It  will  be  noted  that  the  above  equations  contain  the  variables 
)(.,   and    a   from  Fig.   7   which  must  be  expressed  in  terms  of  the 


(13) 
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prlncipal  coordinates. 

In  Fig.  8  by  the  law  of  coain«s 

and  by  the  law  of  sines 

To  sinnplify  Eqs.   (13)  and  (14)  consider  the  following  approximations: 

%,-    Sin    cj^.  Cos    <V,  - 

Then  from  £q.  (13) 

From  £q.   (U) 

^,-5-  %-i':^  [%-%,) 

and  since 

The  approximations  then  are 
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Ji-o      ^Jl 


3  j^^-i^,'^^ 


(15) 


Substituting  <15)  into  (9).   (12).    (11).   and  (10)  yields  the  following  re. 
lations  between  the  extraneous  and  the  principal  coordinates: 


%^  ^^--^ 


il'h) 


(16) 


(17) 


(18) 


(19) 


The  Kinetic  Energy  Term 

The  kinetic  energy  "J"  appearing  in  the  Lagrangian  equation  of 
nv^tion  may  be  derived  as  the  sum  of  kinetic  energies  of  translation 
and  rotation  by  referring  to  Fig.  6: 

Derivation  of  the  Equation  of  Motion  for  the  Coordinate      >^i 

The  first  term  of  the  Lagrangian  equation  may  be  obtained  by 
first  taking  the  derivative  of  £q.  (20)  with  respect  to    ^      : 


il-     M.^^li^      +    M,    i^ 


a<V 


H, 


ill 
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The  extraneous  coordinates    ^-  and    ^    appearing  in  the  above  expres- 
sion niay  be  eliminated  by  differentiating  Eqs.   (16)  and  (18)  and  substi- 
tuting the  results.     Thus  we  obtain: 

The  above  expression  may  be  differentiated  with  respect  to  time  to 
yield: 

The  second  term  of  the  Lagrangian  equation  n^ay  be  obtained  by  differ- 
entiating £vi.   (20)  with  respect  to 

^   =     0  (22) 

The  generalized  force  o(  I    can  be  coxr.puted  by  considering  the 
work  done  by  the  external  forces  as  the  coordinate    Q  ■    is  varied,   all 
other  coordinates  remaining  fixed. 

In  general. 

Workers, n.    ]  °'    Si--   ^^ 

In  particular,   for  coordinate     Q^     ,   the  force  ^      becomes,   referring 

to  Fig.  6: 

^  ^<^^^   '      _       (F  5\r\  c^4  -h  Dcosc^4)  ^  S- ' 

'^        ^%  ~  ^^, 

With  the  approxixnations  of  £q.   (15),   the  above  expression  becomes 

(23) 


5,^     F^^f  D 
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Substituting  £qs.   (21).   (22).   and  (23)  ixito  E^.  (i),  the  equation 
of  motion  in  the     ^.  direction  becomes,   finally: 


Derivation  of  the  Equation  of  Motion  for  the  Coordinate     Vz 
In  a  similar  manner,   differentiating  Eq.  (20): 

and,   differentiating  Eq.  (17): 
Similarly,   differentiating  Eq.  (19) 

Combining  the  above  results: 

By  differentiating  this  expression  with  respect  to  time  the  first  term 
of  the  Lagrangian  equation  becomes: 
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The  second  term  of  the  Lagranglan  equation  is: 

|i     =     0  (26) 

The  generalised  force   S^    ,   referring  to  Fig.   5.   may  be  ex. 
pressed  as: 

With  the  approximations  of  £v^.   (IS) 

5^  ^   ^  (m,^n),)3    -h   F  -D  c^^  (27) 

Substituting  Eqs.   (25).    (26),   and  (27)  into  (1).   the  equation  of  motion 
in  the    ^     direction  becomes: 

(M,.HO^,tKi,r^^,  a --«.)=  F  -3fM,.M.)  -  D<J.^  (28) 

Derivation  of  the  Equation  of  Motion  for  the  Coordinate    >3 

The  first  ternn  of  the  Lagrangian  equation  is  obtained  in  an  iden- 
tical manner  to  that  employed  for  coordinates     ^^   and    ^.  : 
Differentiating  £q.   (20): 

and,   differentiating  £q.  (16): 


I 


h\  M.-^-m,. 
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alao,   by  differentiating  II4.  (18): 


^^3 

Substituting  into  (29) 

u 


L   —     [t    f  H.MxC 
3  p        M.+M,.   _ 


^ 


_       M. Mi5.fr 


M,tM- 


^. 


Differentiating  the  above  expression  with  respect  to  time,   the  first 
term  of  the  L*agrangian  equation  beconnes: 


From  Eq.   (20) 


3     H.  fc^V" 


Q      ^% 


12^ 
0  a^ 


By  differentiating  Eq.   (16): 


'di 


H 


X"      _  _ 


MxM.g. 

(M.  +  Ma)^ 


By  differentiating  Eq.   (18): 


(30) 


(31) 


(32) 


(33) 


"i  >'  1    .      // 

The  term   I  ^-^      which  accounts  for  the  propellent  mass  transfer  out 

of  the  rocket  ia  not  correct,    inasmuch  as  this  moment  should  be  calcu* 

lated  from  the  velocity  with  which  the  propellent  leaves  the  transfer 

tubes  and  enters  the  motor.     Corresponding  errors  appear  in  subse. 

quent  equations.     The  numerical  effect  of  the  resulting  terms  is  negli. 

gible. 
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In  order  to  obtain  ""fe        and       -~- —      ,  It  becomes  neces- 

■ary  to  differentiate  the  exact  expresaiona  for     ^     and    <^     ,    aince  the 

^   term  haa  vanished  in  the  linearisation  of  the  exact  expreaaiona. 
E^a.  (10)  and  (12)  will  be  differentiated  and  the  reaulta  then  linearised. 

From  £q.   (12) 


(12) 


.  }^^  l.-^!  Co 

(34) 


_  M.M^^^4   ^^^^    ,->^^f^^_C0be     ^J^;^cose 


From  Eq.  (13): 
By  differentiation: 


(35) 


7^     - '-^^-^ (2^.^    Or^PrTerrv,; 

^     ~      _     Aj^   r   ct.    -  q  ,  >i  (37) 


f  (%-i3) 
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From  £q.   (14): 

^^   %-  5'"''  li^  ^ '"(%-%)] 

By  differentiation: 

e  =  %  -  -^ 

e  -   i.  -  -It  (%-S)  w 


-f-v 


Alio: 


-^^—    _    cos  6^^^    _    ^^ 


(41) 


(42) 


By  substituting  Eqs.   (35)  through  (42)  into  (34)  and  ignoring  terms  of 
second  order  or  higher  in  the  small  quantities: 

Ak  rr    -  -^^  A  i,    t  ,-^^.  I  %i  (43) 

In  a  similar  manner,  from  £q.   (10): 
By  differentiating  this  expression: 
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sine  bMa. 


Cos  o-  Siwe^^-^ 
^3 


(44) 


^J^llh    C036    -SineP.ll 

Again  substituting  £qs.  (35)  through  (42): 

s:^  -  r^.^M.  '>3  '^(,.,,,M,y  '^^ 

Having  avaluated       _^<l       and       ^^         from  the  exact  expretsiona, 

XT  ^^^  ^'^^ 

^  can  now  be  obtained  by  substituting  these  values  in  £q.   (31). 


(45) 


^^1 


^T 


MiM,^ 


The  generalised  force 
Then,   referring  to  Fig.  6 


(46) 


3  ^^, 
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Substituting  in  values  of    ^     and    ^     from  £qs.   (11)  and  (10) 
Performing  the  differentiation  indicated: 

M.    r    ^ s/H© 


^3  =  [F^'n. -o="v]fe> 


V3 


+  s,-e  J^*  "• 


■ViJ 


t[Fc.^,-Ds.,,][-^,^^.4^.-a-^^. 


^^ 


New,   (ubstltutiag  Eqi.  (3S).   (41),  and  (42): 


^•V; 


1 


+ 


Fiz^ly.    retaining  only  first  order  terms,  the  expression  reduces  to: 

M ,  i, 


^3=^   -^[f(^.-^3)-d] 


(47) 


The  equation  of  motion  for  the  coordinate   Q^     can  then  be  ob. 
tained  by  substituting  £qs.  (30).   (46).   and  (47)  into  £q.  (1). 


IS 


M.  Mii^.  /r  .J       r  -    \     .     mI  K.X, 


M,  ^-Mv 


^^•v^40-7^^BiS^(H^^^iO 
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Derivation  of  the  Equation  of  Motion  for  the  Coordinate     >4 
By  differentiating  £q.   (20)  with  respect  to:  4 

Troxn  Eq.   (16): 

Differentiating: 

From  E4.  (18): 

Differentiating: 

S\abstituting  these  quantities  into  £q.   (4V) 


<49) 


31 


The  first  term  of  the  Lagrange  equation  is  obtained  by  differentiation 
of  the  above  expression: 

From  £q.  (20) 

By  differentiating  £q.   (16): 
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And  by  differentiating  £q.   (18): 


'^7      - 


MiMi-(». 


In  order  to  evaluate      -■      ^  and      -r it  i«  necesaary,   as 

in  the  case  of  the     ^-  coordinate,   to  differentiate  the  exact  expretsion*! 
and  then  make  simplifying  approximations: 
From  £q.   (12): 


cose 


Differentiating: 

By  differentiating  Eq.   (13): 
Also  by  differentiating  £q.  (36): 
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By  differentiating  Eq.   (14): 
And  by  differentiating  E^.  (38): 

le__    Ji.]!UV-'»3)(l-^-iO 

Then       _^ru        nvay  be  evaluated  by  the  substitution  of  the  above  ex- 
pressions  and  £qs.   (35)  through  (42)  into  £q.   (52): 


(53) 


Trom  Ecj.   (10): 

q       r:  q       ■\- -^-^  Jj  .   t^\^  q  q  -  JUL   cob6 --r  '  -  f^ c;^4C05 


Then        •     ■  ^        may  be  evaluated  in  a  similv  manner  to  '^^       by 

substituting  for  the  extraneous  terms  in  £q.  (54): 


The  second  term  of  the  Lagrangian  equation  may  be  evaluated 
by  substituting  the  derivatives  obtained  in  the  foregoing  steps  into 
£q.  (51): 
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The  generalized  force 

Wor  k   <^ 

By  reference  to  Fig.  6,   it  is  evident  that: 

Substituting  for  the  extraneous  coordinates  and  simplifying,   as  in  the 
derivation  of  £q.   (47): 


24-    -t^f^.^^O 


(57) 


Introducing  £qs.   (50),   (51),   and  (57)  into  Eq.  (1),  the  equation  of  mo- 
tion  for  coordinate    ^.  becomes,  finally: 


(Mi 

-_-^  (^.-^0  %  --  '  ^  [-^^.  ^  ^-;  (58) 


-31. 

IV.    COMPUTATION  OF  SYSTEM  COMPONENT 
TRANSFER  FUNCTIONS 

Only  the  two  equation*  involving  the  rotation  of  the  rocket  (co* 
ordinate     ^7  )  *n<i  the  motion  of  the  nozzle  (coordinato     ^  .  )  are  con< 
aidered  since  the  stability  investigation  is  znade  for  rotation  only. 
First,  the  thrust.  Jet  damping,   and  nu>inent  of  inertia  ternns  are  sub« 
stltuted  in  terms  of  the  equation  variables  and  constants. 

Rocket  thrust: 

f  -      _    M,  3    I  sp 

Jet  daxrxping  force: 

(The  negative  sign  is  necessary  in  order  to  make  the  term 
positive,    since    Mi    i*  z  negative  quantity.) 
Moment  of  Inertia  derivative: 

J  I    :::    M»  Jk,^  (^   r  radius  of  gyration) 

•  "2.  • 

I,  -   M.^.  i-2j^,M.-*^, 

X,C:     M.'<t,  (since  ^1    will  bo  small  for 

most  configurations) 

Substituting  the  above  terms  in  Eq.   (48): 


-'^^3h(^i\)-^'^^^,]        (59) 
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Again,   substituting  for  p  ,     D    ,  and    J,    ,  and  introducing  the  applied 
moment  ^  into  Eq.  (58): 


"^J^'^^'^^=  '^^ -f^  i^  fr -f=^.  ^0  -'^ 


(60) 


The  Laplace  Transformation  is  applied  to  the  equations  prior  to 
the  stability  analysis.     Ref.   (2)  explains  in  detail  the  use  of  this  trans- 
formation in  such  analyses.     In  essence,   the  transformation  converts 
the  differential  equations  to  simple  algebraic  equations. 

Since  the  vertical  and  horizontal  velocities    Q,     and    ^     are 
quite  small  during  the  launch,   the  terms  Involving   <^      and   Q    <^      are 
neglected  in  the  equations . 

Application  of  the  Laplace  Transformation  to  £q.   (59)  yields: 


-    Q. 


(61) 


where  Q.    represents  the  transformed  coordinate     Q  .   .   and    S    ^>  the 
complex  variable.     Similarly,   transfornaing  Eq.   (60): 

The  effect  of  the  control  n^oment  iTt  on  the  main  rocket  is  neglected 
in  £q.  (59).  For  a  system  in  which  the  nozale  is  of  comparable  size 
and  weight  to  the  main  rocket,   the  term  should  be  included. 
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^^[L^-^   k.mJ^    ^[(M.tnx)-  ^^  ^'W 


(62) 


where  ^     represents  the  transformed  moment  ^^  . 

In  order  to  solve  for  the  transfer  fiinctione    Fi^    and    Fj    it  is 
necessary  to  rearrange  the  equations  slightly;  Eq.  (61)  becomes 

Similarly,   £q.  (62)  becomes: 


Th«n  the  transfer  function   pL    from  Eq.  (63)  becomes 


F.- 


O3         _    ["mTT?^  ^        k'^^'^O^        M.tM^J-"     LM.tMi        (^M.tMO^J 


(64) 


(65) 


By  substituting  £q.   (63)  into  (64),   the  fvinction    p^,    nrv^y  be  written  as 
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h^ 


>-r <     — 


At  an  aid  in  the  numerical  computations  for  the  stability  analy- 
sis, the  transfer  functions  are  symbolised  as  follows: 

05    _   X 


F,= 


<?4-Q-i 


Y 


I 


M 


a;2,-w 


(67) 


(68) 


where: 


_  r[i:},!:kMl  ,v  [wv^Mj.K   M>ai£)     Unigirp   MXii^l] 


Wr 
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V.     SYNTHESIS  OF  THE  NOZZi-E  CONTROL  TR/\NSFER 
FUNCTION   "Fi".   .\ND  STABILITY  INVESTIGATION 

To  illustrate  the  determination  of  the  nozsle  control  transfer 
function,   a  specific  numerical  example  will  be  worked  out  using  the 
German  V-2  as  a  typical  vehicle. 

As  indicated  in  the  Analysis,   the  system  transfer  function  may 
be  written  as 

Also,   the  computer  synthesizes  the  dynamics  of  the  rocket  to  give 

/^  =  >>  (70) 

Substituting  Eq.  (70)  into  (69) 

^.6.=        .  ^^/"  ^i (71) 

Let 


(72) 


Substituting  for    /^    and  /^     from  Eqs.   (67)  and  (68): 

rs  =/-  / — ^ — ; 

^  (  yz  -  y^/  J 

AS  a  numerical  example,   consider  the  V-2  rocket  with  pivoted 
rocket  motor  and  nozzle.     The  paramieters    X    >     Y    »    Z    >   ^^^  ^ 
have  the  following  values  as  calculated  in  the  .ippendix: 


(73) 
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y=   92,  3 ao  s"-  -    no  5 

2=       Zo  3  O    3^     -    //.  (>S  S 
^=      6'^-^    S'-      V-       /o    ff  S  S 


> 


(74) 


The  problem  now  is  to  determine  a  function    /-^    which  will  give 
the  systenrn  suitable  stability  characteristics  during  the  launch.     The 
procedure  will  be  to: 

1.  Assume  a  function    f-,     . 

2.  Check  the  system  for  stability  by  the  Nyquist  criterion. 

3.  Having  established  a  stable  systen^,   determine  the  roots  of 
the  system  transfer  function  by  means  of  the  Evans  Root 
Locus  Method. 

4.  Knowing  the  roots,   at  a  typical  operating  gain  of  the  feedback 
system,  investigate  the  damping  of  the  system. 

5.  If  the  resulting  motion  is  unsatisfactory,   improve  the  choice 
of    /^    and  repeat  the  process. 

In  explanation  of  the  methods  to  be  used  to  investigate  the  sta- 
bility of  the  feedback  control  system,   it  will  be  noted  that  the  equations 
of  motion  are  linear  differential  equations  with  constant  coefficients. 
The  solution  for  the  particular  coordinate  in  question  can  then  be  rep- 
resented in  the  general  form 


%-Q^   +c 


2^       ^ 
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wh«re    5      is  a    siziiple    root  of  the  •quation.    In  order  to  have  a  damped 
motion  which  will  subside  with  time,  it  is  obvious  that  the  complex 
roots    s-     must  all  have  negative  real  parts,      uiother  way  of  stating 
this  requirenient  of  a  stable  system  is  that  all  the  roots  s^      of  the  sys- 
tem must  appear  in  the  left  or  negative  half  of  the  complex  plane. 
The  Nyquist  Criterion 

The  Nyquist  Criterion  is  an  application  to  the  simple  feedback 
network  of  the  Cauchy  Theorem  in  the  theory  of  complex  variables. 
This  theorem  cnay  be  employed  to  show  that  if    6/51     is  the  vector  rep- 
resenting the  function  of  a  complex  variable  ^5    ,  the  function  having 
/77    poles  and  n    aeros  within  a  closed  contour  C    in  the    s    plane,  that 
as  the  point    S    moves  around  the  contour  once  in  a  clockwise  direction, 

the  vector     o^^.     carries  out    n-m    clockwise  revolutions  about  the  ori- 

//   i> 
gin.    If  the  contour    C    is  chosen  to  include  the  entire  right  half  plane, 

and  the  function  in  question  is  plotted  as  the  complex  variable  encircles 

the  right  half  plane,  the  difference  between  the  number  of  seros  and 

poles  of  the  function  that  exist  in  the  right  half  plane  will  be  indicated 

by  the  nunriber  of  encirclements  of  the  origin.     Either  the  number  of 

seros,  or  poles  existing  within  the  contour  must  be  known  previously. 

In  applying  the  theorem  to  the  feedback  network,  the  function  in 

question  is  usually  expressed  as 

where  the  function      n^(S         can  be  factored;  thereby  determining  the 
zeros  of  the  function    n^^^  The  denominator    /-^  ^G     in  general 

cannot  be  factored  so  readily,  but  application  of  the  Cauchy  Theorem 
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will  show  the  number  of  roots  of  the  denozx\izM.tor  existing  in  the  con- 
tour chosen,   or  the  right  half  plane. 

In  order  to  obtain  a  plot  independent  of  the  gain      ^ ,   the  inverse 
of  the  function  may  be  considered. 


/     _ 


-    Z^    -/ 


/Tg     65  /f    6 


/ 


/ 


Then  if  the  vector      ^    g         encircles  the  origin,   the  vector    ^p^ 
will  encircle  the  point     -/  ,   as  the  variable    s   moves  around  the  con- 
tour of  the  right  half  plane.     Or  sixnilarly  the  vector    ■q'       will  en- 
circle  the  point       ^  .     The  function      -q~      is  then  plotted  as     s  en- 
circles the  right  half  plane.     By  application  of  this  Nyquist  Criterion, 
the  stability  or  instability  of  the  system  can  readily  be  determined. 

If  the  noszle  control  transfer  function   ^^     is  first  assumed  to  be 
equal  to  a  constant  f '^  ~   ^/)     ^^  system  transfer  function  takes  the 
form 


*    '"    7775"    °      ...r^JL.     -,  (T5) 


For  the  Nyquist  diagram,   the  inverse  of    G    is  plotted  as  the  vec< 
tor   5    encircles  the  right  half  of  the   S    plane.     The  function  to  be 
plotted  then  is 

—^  =  —X —  <"» 

In  order  to  use  the  Nyquist  stability  criterion  it  is  necessary  to  know 
whether  the  function    ^G     has  any  seros  in  the  right  half  plane.     By 
factoring  £)q.   (73)  the  following  expression  is  obtained: 
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^G^      =    ^jl-nix/O  ^   )     ^ ^ ^ (77) 

Since  the  function  V^  "h*«  one  »ero  In  the  positive  real  hiOf  of  the 

_/ 
plane«   the  Nyquist  Criterion  requires  that  the  plot  of    q    '  must 

encircle  the  point  -4^  once  in  a  counter-clockwise  direction  as  the 

S   vector  circles  the  right  h^lf  of  the   s    plane  once  in  a  clockwise 

direction. 

Fig.   9  is  the  Nyquist  diagram  of  Eq.   (76).     As  the  frequency 
increases  from  sero  to  positive  infinity,  thence  to  negative  infinity  and 
back  to  sero,   there  is  one  clockwise  encirclement  of  all  gains  from 
^  to  oo   .     The  net  encirclement  may  be  considered  to  be  two  clock- 
wise rotations,   indicating  two  roots  of  the  system  transfer  function 
existing  in  the  right  half  plane.     The  system  is  unstable  for  all  gains. 

Next,   assume  the  nozale  control  transfer  function  to  be  of  the 
form  f^  =  K,  0^  7S)  ;  then 

/Cs  Os   --       ii^J^lM^A.^ (78) 

The  Nyquist  diagram  for  the  above  equation  is  shown  in  Figs.   10  and 
11.     Again  the  number  of  eeros  of    ^6    existing  in  the  right  half  plane 
determines  the  Nyquist  criterion.     The  factored  equation  now  becomes 

^^rsj'-^M^^^^K? — -^, ^^^ >f  <^^> 

Since  there  is  one  zero  in  the  right  half  plane,   again  it  is  desired  to 


//,   /« 


have  one  counter-clockwise  encirclem.ent  of  the  gain  o     as  an  indica- 
tion of  sUbility. 
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It  was  found  necessary  to  plot  the  Nyquist  diagram  to  several 
scales  to  indicate  clearly  the  crossing  points  of  the  negative  real  axis. 
The  enlarged  portions  were  plotted  for  the  case  when    7=/       ,   and  ap- 
pear as  Figs.    12  and  13.     From  these  plots  it  is  clear  that  there  is 
one  coxinter -clockwise  encirclement  of  ail  gains  fronn    A^^  /^  to 

^-  75"  .     The  configuration  is  stable  for  gains  in  this  range. 
The  Evans  Root  Locus  Method 

This  niethod  is  a  graphical  solution  for  the  variation  of  the  roots 
of  the  feedback  system  transfer  fvmction    /Y  G       ^^  the  gain   K  is 
varied  from  zero  to  infinity,  where 

The  plot  is  constructed  by  observing  the  following  rules  (Ref .   2): 

1.  For    nz:Q     ,   the  poles  o(    /fl  ^     aire  the  poles  of  (^ 

2.  For  4<-^m  c^  ,  the  poles  of  /f"   q*    are  the  aeros  oi  (^  , 
with  any  missing  poles  taken  at    s  -  a:^  » 

3.  The  product  of  the  distances  from  a  root  to  the  poles  of  /f^ 
divided  by  the  product  of  the  distances  from  the  root    to  the 
seros  of    A'G    is  equal  to  the  gain    4^  I' 

4.  The  sum  of  the  angles  from  the  aeros  of    /TC    to  a  root,  less 
the  sum  of  the  angles  from  the  poles  of    ^^    to  the  root  must 
equal  7?    ,  or   7?  plus  an  integer  nnultiplied  by     ^  Z?'. 

The  locus  of  the  roots  of  /f^  ^  can  be  determined  graphically 
by  trial  and  error,  observing  the  above  rules.  The  gain  scale  can  be 
plotted  along  the  path  of  the  roots,   and  for  a  chosen  operating  gain. 
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the  roots  can  be  determined  by  inspection. 

The  root  locus  plot  of  Eq.   (18)  for  the  case  when    7=^  /      is  shown 
in  Fig.    14.     The  roots  emerge  from  the  poles  of   A^G    at  aero  gain  and 
proceed  along  the  loci  indicated  to  the  seros  of   ^ C    at  a  gain  of  infin- 
ity.   Gains  at  intermediate  points  can  be  approximated  by  dividing  the 
product  of  the  distances  to  the  p>oles  of   ^G    by  the  product  of  the  dis- 
tances  to  the  seros,   and  by  the*constant  gain  factor  in  the  equation 
(49.2  X  10~   ).     It  is  apparent  fron\  the  root  locus  plot  that  the  range 
for  a  stable  configuration  is  from  the  gain  where  the  roots    --^  3 
cross  into  the  stable  half  plane  to  the  gain  where  the  roots  5^  ^ 
cross  into  the  unstable  half  plane,   referring  to  Fig.   14.     The  loci  in 
this  range  are  plotted  to  an  enlarged  scale  and  appear  in  Figs.   15  and 
16.     As  an  example  of  the  computation,   approximate  expressions  for 
the  gain  along  the  loci  in  this  range  are: 

For  root 

^      _,    r'^^ZS'-  W<^'^  ^^<^J  ^coV ^  ^  ,.    .. 

For  root 


rJ/J97^/)f/s?.  usps-yi^sf.  /^J9/)^9.c  ^^^"^ 

It  will  be  noted  that  the  gains  at  the  limits  of  the  stable  range 
check  within  limits  of  accuracy  of  the  scale  used  on  the  graph  as  shown 
in  the  table  below  for  the  Ny^uist  and  the  Evans  Root  Locus  plots. 
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Method 

Upper  Limit  of  K 

Lower  Limit  of  K 

Nyquitt  plot 
Evans  Root  Locus 

75.0 
74.80 

12.0 
11.85 

Having  established  the  gain  scale  along  the  root  loci  in  the  range 
of  positive  stability,   referring  to  Figs.   15  and  16,   it  is  possible  to 
choose  an  operating  gain  and  determine  the  values  of  the  roots   ^3 
and  S]^  y      at  that  gain.    Some  possible  operating  points  have  been 
tabulated  below 


Roots 

K  3  50 

K  >  60 

K  •  70 

Si 

-.0031  +  i  .156 

-.00325  +  i  .171 

-.0033  +  i  .185 

S3 

-.0031  - i  .156 

-.00325  -  i  .171 

-.0033  -  i  .185 

S2 

-.0068  4-  i  6.25 

-.004  -f  i  6.25 

-.0003  ^i  6.25 

S4 

-.0068  -  i  6.25 

-.004  -  I  6.25 

-.0003  -  i  6.25 

The  low  frequency  roots  Sy ^        are  associated  with  the  motion 
of  the  main  rocket,  whereas  the  high  frequency  roots    S^  y     are  as- 
sociated with  the  motion  of  the  nossle.     The  time  to  damp  to  half  am- 
plitude and  the  period  may  be  computed  for  the  modes  quite  simply; 
the  results  of  such  a  computation  for  a  gain  of    ^=^0  are  tabulated 
below. 
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K  »  60 

Rocket  Oscillation 

Nossle  Oscillation 

Time  to  damp  to 
half  amplitude 

210  seconds 

171  seconds 

Period  of  one 
cycle 

36.8  seconds 

1  second 

Number  of  cycles 
to  half  amplitude 

5.7  cycles 

171  cycles 

It  is  apparent  that  the  damping  is  quite  poor,   and  an  improve- 
ment is  desirable.     The  acceleration  of  the   V-2  at  take-off  is  roughly 
equal  to  30  ft.  /sec.^.    If  it  is  assumed  that  at  a  velocity  of  200  ft.  /sec 
the  aerodynamic  forces  stabilize  the  rocket,  the  period  of  take-off  in- 
stability is  approximately  7  seconds  long.     For  the  foregoing  exanr\ple, 
the  rocket  conceivably  could  complete  only  one -fifth  of  a  cycle  of 
slightly  damiped  motion  in  the  take-off  period  as  a  result  of  a  unit  dis- 
turbance at  the  instant  of  firing.     The  noaale  oscillations  are  so  rapid 
as  to  conceivably  cause  overloading  of  the  control  device. 

In  order  to  improve  the  damping  characteristics  of  tho  system, 
referring  to  the  root  locus  plot.   Fig.   14,  the  addition  of  one  or  more 
aeros  along  the  negative  real  axis  would  serve  to  pull  the  roots    $, 
and    S-T  into  the  negative  half  plane  more  quickly,  since  zeros  attract 
the  roots  much  as  lines  of  force  are  drawn  to  nr\agnetic  poles.     This 
should  increase  the  magnitude  of  the  negative  real  part  of  the  roots  at 
a  given  gain.    vUso,  the  angle  of  departure  of  the  root    Sp    would  be 
increased  counter-clockwise  as  a  resxilt  of  the  angle  condition  of  the 
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root  locus  plot.     The  root  would  then  be  forced  to  travel  a  longer  path 
before  passing  into  the  unstable  half  plane. 

The  addition  of  seros  in  the  left  half  plane  involves  adding  addi* 
tional  lead  networks,   which  in  real  connponents  have  transfer  functions 
of  the  type     (/^^js)      where     (f^^-^Ti^ 

Assuming  a  new  nozale  control  transfer  function  to  be  of  the  form 


f'    >» 


where     7~-  /     ,   the  factored  expression  for    /CQ    i 


s  now 


(83) 


The  root  locus  plot  of  £q.  (33)  is  shown  in  Fig.   17.     The  range 
of  gain  for  stability  is  now  from  Ks4toKB913. 

By  approximating  the  loci  of  the  roots  in  the  stable  range  by 
straight  lines,   the  roots  associated  with  various  gains  can  be  estimated 
conservatively.     Some  representative  values  are  tabulated  below. 


Root 

K  =  60C 

Si 

-.0446  -f  i  .541 

"^3 

-.0446  -  i  .541 

Sz 

-.02  +  i  15.25 

24 

-.02  -  I  15.25 

If  the  damping  time  and  period  are  computed  for  the  latter 
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•ystexn,   the  following  results  are  obtained. 


K  •  600 

Rocket  Oscillation 

nozzle  Oscillation 

Time  to  damp  to 
half  amplitude 

15.3  seconds 

34.2  seconds 

Period  of  one 
cycle 

11.6  seconds 

.412  seconds 

Number  of  cycles 
to  half  amplitude 

1.3  cycles 

83  cycles 

Interpreting  the  results  as  before,   the  rocket  now  completes  ap- 
proximately two -thirds  of  a  cycle  of  oscillation  during  thm  launching 
period,  but  the  effect  of  a  unit  disturbance  is  damped  to  a  much  greater 
extent.     The  nossle  oscillations  occur  at  a  higher  frequency  but  the 
damping  has  been  improved  considerably. 
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VI.     CONCLUDING  REMARKS 

The  Mialysis  shows  the  feasability  of  stabilizing  a  rocket  during 
the  launching  period  by  utilising  the  motion  of  a  compound  pendulum 
suspended  nossle  as  the  input  to  a  feedback  servo. control  system.    It 
was  determined  that  the  nozzle  control  for  such  a  system  must  include 
response  to  both  the  attitude  and  the  rate  of  change  of  attitude  of  the 
rocket,   as  computed  from  the  pendulum  nvotion,   in  order  to  obtain 
stability.     It  was  further  determined  that  changing  the  nozzle  control 
to  include  response  to  accelerative  motions  of  the  rocket  improved  the 
daxnping  characteristics  of  the  system. 
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APPENDIX 


The   V-Z  Rocket 


Soo    Kg 

-1- 

1 

1 
1 
I 
1 

/ 

/ 

-   + 

3SOO    Kj            1 

f  e  ; 

5250  K^ 

1  ^ 

i  4 

— 7  42'h. 

AUKILIARieS       1 

1 
■ 1 

ALCOHOL        [ 

OAV*rN 

' PUMPS 

1 

'            MOTOR 

\ 

20  83 

—         37 

s' 

JO. 2 

^/   7        

OVERALL       HeiCHT 

Weight  and  Balance 

jo' 
Table 

— 

Wt.  in 

K£.. 

Wt.  in  Lbs. 

Moment  Arm,  ft. 

Moment 
lb.  ft. 

1.  Explosive 

980 

2160 

7.42 

16,000 

2.     Fuselage 

1750 

3860 

26.65 

103,000 

3.  Pump 

450 

993 

37.5 

37,200 

4.   Motor 

550 

1214 

41.7 

50,  600 

5.  AuxilUrie 

s        300 

662 

14.0 

9,300 

6.  Alcohol 

3500 

7720 

20.83 

161,000 

7.  Oxygen 

5250 

11600 

30.2 

250,000 

Totals 


12,980 


28,  620 


727, 100 


Moment  arm  of  C.G.   =  ^^^*^^^  =  25.4  ft. 

28,620 


Note:    Data  was  obtained  from  Ref .  (4)  and  is  onl^  approximately 
correct. 
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C&lculation  of  the  Moment  of  Inertia 


8  75 


h'3S  35 


4-    CG 


25  4i 


\ I 


Aituxne  the  rocket  is  a  solid  cylinder  at  shown  in  the  above  sketch. 
Radius  of  gyration  of  a  solid  cylinder  about  its   C.C.  is 

>3  --    V  — I 


'OK 


^:^^"-^^ :  vW  r  ^.7^ 


4t 


Additional  DaU  from  Ref.   (4) 


F  s  Total  Thrust  =  27.  2  tons  ^  54.  400  lbs. 


W  s  Total  Weight  »  12.  980  Kg.  »  28.  620  lbs. 


Mj  .  Mass  Flow  =-  iiii£2i^  =  _  193.  2  lb. /sec.  .-6   lb. sec. /ft. 

100  sec. 


l-o  «  Specific  Impulse  =  -r-^  =  51|J00   ^  ^82  sec. 
■P  Ml  g       193.2 
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Calc Illations  for  the  Constants  of  S«ction  III 


M.  =    Z7.  406       lbs     -     851 


lb  sec 


ft 


M,  --    I,ZI4    Ib5    =      37.  y  -^^^ 

k,  -    9  V    ff 

[^^  o^      2.      ft       (rA<imS       of     jyritiorx      of  irtotor) 

;,  =  2-0  f^ 

i,.  -    3   7    ^t 

A  '    0.3    ft 

is  -  4    0   ft 

I,  =    M,  k,^     -     80^600      lb    ft    s<?e^ 


I^  -    Mtk^""  -     /5I      lb  ft    sec"^ 


1 


PIN 


X  = 


M.M.i.i.  z    .        /M^2lAA      ^    H.M.X./.\  /n.M,/gl5r        ^V^'>^.A] 


X- 
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(Z7^^06  lks)(i77—iv^)  do  fi)(^.7ff) 


^I  6^o  \bs 


S'     -H 


(83  8   7    ib.Z££Mi  Bftft   7    [^-ISs} 

^  Pr       '  fir 


(s^r-HP  )(-(So'^|^)(^-W(5^.2U.)^8^  -c; 


9ftft.7 


^     ^      [26  75     s'  -    4<;/     S    -      I  o<^^     X    10*  1        (»6/0 


Y=    (I, 


n.M^^.(i'.-i',)\  ^2. 


M,  ^f-l^ 


Y  = 


■(ao,aoo  „,..,., /e.i^)^3..^^^;.)a.-3.;^>^  ^. 


(-.0  v^)(.7../ .  (n:z^'^^(-^-^^y^'^^)(^o.,,,,. 


( 


(8^'  ^  Ve.o  ^^  )(zofi)Uo(f) 


86B.7     tk^S5S-^ 


Y   =    [  9Z,  3ao  s*-  -  //os] 


Z  = 


z- 


fe5l)(J7  7)(3  ?)(20-3  7j 
880  7 
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I5W  5^  f     ^^^^^^^"^^^^^^^^^~^  ''^ 


(8S  8  7;*^ 


^«'<^)(-^^^^;; 


7    =    [  Z030   s'   —    )».6S     5    ] 


\A/  = 


(x.^^^).'-  (^.  -"•^(^■-7^/.))^- 


1  •  »■ 


w  = 


(-  ^'^'^'^^)^'-(^'^l^^^-i-^H<^-'^',^))s 


(:i7  7)\-6)^ (10 )(i  7) 


(868  7j* 


\a^    =      (iy»   -*-495^  5*"     -*-     (-/5  -h  n )    5     —.005-4- 


\a/    ^      [646   S^     -».      10  85"    S  J 
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